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We present a new approach to analyze homogeneous nucleation based on non-equilibrium ther- 
modynamics. The starting point is the formulation of a Gibbs equation for the variations of the 
entropy of the system, whose state is characterized by an internal coordinate or degree of freedom. 
By applying the method of non-equilibrium thermodynamics we then obtain the entropy produc- 
tion corresponding to a diffusion process in the internal space. The linear laws together with the 
continuity equation lead to a kinetic equation of the Fokker-Planck type. By choosing properly the 
degree of freedom we are able to obtain a new kinetic equation for a global crystallization order 
parameter ( used in recent simulations ), and also we recover some of the existing equations. The 
consistency of the scheme we propose is proved in the quasi-stationary case. Finally, we also outline 
<D . the way in which our formalism could be extended to more general situations. 
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> . I. INTRODUCTION 

on 

The first step in many phase transformations of pure substances -for example condensation of gases or solidification 
j of melts- is the formation of small embryos of the new phase within the bulk metastable substance. This is an 
activated process: a free energy barrier must be overcome in order to form embryos of a critical size, beyond which 
the new phase grows spontaneously. This fundamental mechanism of phase transformation is known as homogeneous 
ON . nucleation, and the rate at which critical-sized embryos are formed is the nucleation rate ( for reviews on nucleation, 
' see for example references |jj - || ). 

To analyze the process, three different types of schemes have basically been proposed in the literature. The first one 
is referred to as classic nucleation theory. In it, the process is usually described by means of a master equation that 
completely determines the evolution of the clusters size distribution. The fact that the size of the clusters vary in a 
discrete way (due to the lost or gain of a single monomer) causes difficulties in the treatment and resolution of this type 
of equation. For this reason there have been several trials to obtain a continuous diffusive equation as an approach to 
the master equation. The advantages that a continuous description may offer are mainly the computational efficiency 
and the possibility of transforming the Fokker-Planck equation into a Shrodinger-type equation, soluble by means 
• i-H , of Quantum Mechanics' usual methods. The second type of description is centered on the construction of dynamic 
models of the cluster interface. By assuming that the molecules of the interface behave as Brownian particles, it 
is possible to arrive at the characterization of the clusters size distribution in terms of a Fokker-Planck equation. 
Finally, the last scheme is based on the theory of the density functional. In this approach to the problem, the kinetic 
description is carried out through the continuity equation for the probability density. 

As we will see, the three points of view lead to express the nucleation rate by means of a diffusive equation. However, 
in spite of the great variety of proposed equations, up to now none of them has turned out to be fully satisfactory. 

In this paper we present an alternative approach that allows us to obtain kinetic equations applicable to homogeneous 
nucleation. Our method is based on the formalism of internal degrees of freedom Q- M, formulated in the framework 
of non-equilibrium thermodynamics, and applied successfully before to other problems like the adsorption kinetics 
of particles on surfaces jl O fl , the description of chemical reactions |ll]] or the treatment of noise in interfaces of 
semiconductor devices |12|| . This formalism allows us to derive Fokker-Planck equations from a Gibbs equation and 
not as a continuous approach to a master equation. Moreover, within this framework we recover directly the proper 
equilibrium distribution and we may derive explicit expressions for drift and diffusion coefficients. 

The paper is structured as follows. In section II we will discuss the different procedures proposed (inside the scope of 
homogeneous nucleation) to obtain Fokker-Planck equations. Section III is devoted to present our method to derive a 
general Fokker-Planck equation describing the kinetics of homogeneous nucleation. In section IV, we will particularize 
for several different interpretations of the internal coordinate. We will derive a new kinetic equation in terms of an 
order parameter associated to the degree of crystallization of the system and we will reproduce some of the kinetic 
equations proposed in the literature. In section V, we will discuss the quasi-stationary case and show how in this case 
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our approach overcomes some of the inherent difficulties of previous treatments. Finally, in the last section we will 
stress our main results. 



II. BRIEF REVIEW OF PREVIOUS TREATMENTS 



In this section we will present a brief summary of the different Fokker-Planck equations that have been proposed 
to treat the problem of homogeneous nucleation. These equations have been obtained mainly from three different 
ways: as continuous approaches to a master equation; starting from microscopic kinetic models by considering the 
separation of a molecule from the cluster surface like a problem of Brownian diffusion through a barrier; or as the 
consequence of imposing a continuity equation for the probability density in the framework of the density functional 
theory. We will discuss all of these in different subsections. 



A. Continuous approaches to a master equation 

In the classical model of nucleation || one assumes that the size of a cluster can vary due to the lost or absorption of 
a single monomer. Hence the evolution of the clusters distribution is described by the set of discrete master equations 

= p( n -l)p(n-l,t) + a(n+l)p(n+l,t) - p(n,t) (o(n) + (3(n)) (1) 

at 

where p(n,t) represents the density of clusters containing n molecules and a(n) and P(n) are the rates at which a 
cluster of size n loses or gains a single molecule, respectively. 

To transform this master equation into a continuous equation, several different methods have been proposed in the 
literature. The most straightforward one consists of replacing the discrete variable n by a continuous variable x. After 
expanding in Taylor series we arrive at the well-known Kramers-Moyal expansion, which upon truncation in second 
order yields the diffusion-like equation 

% = I (dw£ + W) (2) 



dt dx \ dx 

where A(x) and D(x) are the drift and diffusion coefficients, respectively. 

This equation, however, has the serious inconvenience that it does not reproduce the equilibrium distribution of the 
original master equation. Additionally, it has also been proven that it overestimates the nucleation rate exponentially, 
with an error that grows unboundedly for large n fll3f . 

The Fokker-Planck equations of Frenkel-Zeldovich Goodrich [51| and Shizgal-Barret Jll| try to solve these 
difficulties by imposing a relationship between drift and diffusion coefficients in order to ensure the correct equilibrium 
distribution. These three proposals differ in the expressions for those coefficients, chosen to improve the approximation. 
However, as pointed out by Gitterman and Weiss |l7j ], none of them seem to be able to reproduce the correct 
deterministic behavior . 

A variant of these equations is the one formulated by Rabin and Gitterman |l§| ], where the innovation consist in 
considering that the coefficients of the Fokker-Planck equation cannot be fixed arbitrarily, but rather - in the vicinities 
of the critical point - they are determined by the critical dynamics of the system. 

Van Kampen |19|j proposed an alternative method, consisting in rescaling the rate constants a and j3 in terms of 
a parameter f2 (usually the system size ) and to expand the result in powers of In this way one obtains an 

expansion similar to Kramers-Moyal that shares the same problems described previously. 

The contribution of Grabert, Hanggi and Oppenheim and later Hanggi et al. |Hj is also based on a rescaling of 
the rate constants, but they obtain the coefficients (expressed in terms of the rates a and (3) from non-linear transport 
theory. Their equation reproduces the equilibrium distribution appropriately, but it does not describe the short time 



behavior correctly and it gives an erroneous value of the variance of the distribution |17|. 

Finally, inside the scope of Fokker-Planck equation obtained starting from the master equation, we can quote the 
essentially mathematical method of Wu Q , p2[ that proposes an index change and a mathematical optimization to 
minimize the error associated to approach the master equation - discrete - by means of a - continuous - diffusive type 
expression. 
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B. Kinetic Models 



A different approach to obtain Fokker-Planck type equations for homogenous nucleation, is based on the peculiar 
dynamics of a monomer near the surface of a cluster. Inside this scope we will consider the works of Lovett p3[ and 
Ruckenstein and Narsimhan and Nowakowski p3] - |^7j . 

Lovett uses Thermodynamics to estimate the radius of the critical cluster and the height of energy barrier that 
a monomer should surmount to leave the cluster surface. Considering that the monomer undergoes a Brownian 
movement, he uses the corresponding Fokker-Planck equation to study the probability of leaving the cluster. His 
model is only applicable to nucleation in gases, and he uses the macroscopic concept of interfacial tension which can 
doubtfully be applied to the description of small clusters. 

Ruckenstein and coworkers (24|] - |2^] propose an alternative approach based on the calculation of the rate a. Their 
starting point is the estimation of the potential barrier that is generated at the surface of a spherical cluster of radius 
R from the interaction between monomers and the recount of the number of monomers that interact with a given 
one. They also consider that the movement of a molecule is Brownian. In the case of nucleation in liquids, they write 
a Smoluchowski equation in spherical coordinates for the distribution of probability p{r , t) 

dp{r,t) / - . , D 



V r - • [ DV FP (r, t) + —p(r, f)V^(r) ) (3) 



dt V rrv ' ' kT 

Here D is the diffusion coefficient and <f>(r) the interaction potential of a surface monomer. For gases, they use a 
diffusion equation in the energy space. These equations are just employed to obtain the rates a and f3. The rates 
are introduced in the common Kramers-Moyal truncated expansion of the master equation to evaluate the nucleation 
rate. 

The main inconveniences of this formalism are that it considers (in the case of liquids) that the diffusion coefficient 
is constant, when in fact there are reasons to think that it is not totally true; moreover the contribution of the 
bulk medium to the effective potential is not incorporated rigorously |Bj . On the other hand, an important difficulty 
linked to this type of descriptions arises from the fact that the nucleation rate is very sensitive to uncertainties in the 
interatomic potential p8[. 



C. Density functional theory 

One of the main drawbacks of the classic theory is the assumption that the interface of the cluster is sharp, 
hypothesis that - at least at the level of simulations - is not completely satisfactory. 

Field models allow a treatment keeping in mind the possibility that the interface between the cluster and the original 
phase may be diffuse. The theory of the density functional [Q , |28| considers the free energy as a functional of the 
density and uses it to built a coarse-grained free energy. The free energy barrier can be taken out from variational 
principles. However, to be able to obtain kinetic information, hydrodynamic theory should be introduced. 

The starting point of Cahn and Hilliard kinetic treatment [ p9| , |30[| is the formulation of an equation of continuity 
for the current of probability density of a certain configuration of the system. Following the scheme of non-equilibrium 
thermodynamics, they postulate the following linear relationship 

OF 

J = -MV%- (4) 
op 

where J is the current, M is the atomic mobility (taken as constant), p is the probability density and F is the 
coarse-grained free energy. With these ingredients a Fokker-Planck equation is obtained. 

The work of Langer |3lJ - |?3| constitutes an extension of the previous treatments in the sense that it provides 
Statistical Mechanical basis where fluctuations are taken into account, in the sense that he considers diffusion in the 
phase space. He proposes the following equation 

d P ({y},t) ^ d ( ( dp dOM) P \\ ... 

where {7/} is the set of all degrees of freedom that characterize a system configuration, G({r/}) is the coarse - grained 
free energy functional, and My is a generalized mobility matrix. 

The main inconveniences of these theories based on the functional density are that their correctness depend on the 
appropriate election of the free energy functional and that its range of application is restricted; concretely they are 
not easily applicable to time dependent situations. 
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III. NON-EQUILIBRIUM THERMODYNAMICS APPROACH. INTERNAL DEGREES OF FREEDOM 



In our approach to the problem of homogeneous nucleation, we consider that the state of a system may be charac- 
terized by an internal coordinate or degree of freedom 7 , which may for example represent the number of monomers 
of a cluster (conveniently normalized to become a continuous variable), its radius, or even an order parameter. The 
nucleation process can be viewed as a diffusion process through a barrier, in the space spanned by the values of the 
internal coordinate, from an initial value 71 , corresponding to the metastable state, to 72 , characterizing the new 
phase. 

Our starting point will be to assume that the variations of the entropy, SS , due to the diffusion process are given 
by the Gibbs equation 

1 P 2 

5S = -- fx( 7 ,t)Sp( 7 ,t)d 7 (6) 

•'71 

Here p(j,t) is the density in the internal space, p,(j,t) its conjugated chemical potential and T the temperature of 
the system which is assumed constant. 

In the absence of interactions among constituents of the system with different 7's, the chemical potential is the one 
for an ideal system |m| 

A*(7,t) = fcsTlnp( 7 ,t) + $(7) (7) 

where $(7) is the potential and ks the Boltzmann constant. 

Quite generally, the evolution of the density in the internal space is governed by the continuity equation 

dp _ dJ 

at ~ (8) 

which introduces the diffusion current in the internal space J (7, t). The expression for the entropy production related 
to the diffusion process can be obtained from ^ by using (||). One then arrives at 



when a partial integration has been performed, with the assumption that the diffusion current vanishes at the initial 
and final states of the process. 

Assuming locality in the internal space, for which only fluxes and forces with the same value of the internal 
coordinate are coupled, one then derives from expression (t|) the linear law 

JM . -M!m a») 

when £(7) is a phenomenological coefficient, which may in general depend on the internal coordinate. This expression 
can be used in the continuity equation (^) thus arriving at the Fokker-Planck equation 

8 ^(7,*)^ +6(7,*)^ (ID 



dt d"f \ 97 

where 



is a mobility in the internal space and 



K7,t) = (12) 



D( 7 ,t) = k B Tb(j,t) (13) 



the corresponding diffusion coefficient. 

In order to fully determine our Fokker-Planck equation we need to specify the form of the diffusion coefficient in 
the internal space and the potential. To obtain a microscopic explicit expression for the phenomenological coefficient 
£(7), we will apply fluctuating hydrodynamics in the internal space, in a similar way as in rcf. |DJ]. Adopting this 
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framework, the presence of fluctuations in the diffusion process along the internal coordinate can be taken into account 
by adding a random contribution J r (j,t) to the systematic contribution J s (7,i) in the diffusion current 



J( 7 ,f) = J s (l,t) + J r (j,t) (14) 
Considering that J r (7,i) has zero mean and satisfies the fluctuation-dissipation theorem 

(J r ( 7 ,t)J r (7',iO) = 2fc fl L( 7 )a(7-7 / )*(t-0 (15) 
one can obtain the following Green-Kubo expression for £(7) 
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L{l) = — b \ H\ *(J r (7,0)J r (7',i)} (16) 



Hence, in view of eqs.([l2|) and ( |l3| ) the diffusion coefficient D(j,t) is given by 

1 r72 poo 

D( 7 ,t) = -^J h j o dt(rfr, o)rtf, t)) (17) 

This equation offers an explicit expression to determine the diffusion coefficient in the internal space. For each 
interpretation of the degree of freedom, the later equation could be transformed in order to be evaluated by means 
of simulations. For instance, the way to calculate the diffusion coefficient when the internal coordinate represents a 
global order parameter will be extendedly explained in a future paper. 

Additionally, the nucleation barrier $(7) can be estimated by means of simulations in the way proposed by van 
Duijneveldt and Frenkel in ref. |3^]. Consequently, our scheme becomes in this way completely soluble at least at the 
simulation level. 

Finally, it is worth pointing out that our formalism recover in a natural and direct way ( without additional 
impositions ) the correct equilibrium distribution. Since at equilibrium, the current in the internal space J e „ vanishes, 
from equation (|l(]), this condition is accomplished when the chemical potential is constant. Using eq. (Q) we then 
obtain the proper equilibrium distribution 

p e «( 7 ) oce-*^/ feBT (18) 



IV. PARTICULAR EXPRESSIONS OF THE FOKKER-PLANCK EQUATIONS. 

The formalism we have proposed enables us to formulate a general Fokker-Planck equation in terms of an arbitrary 
internal coordinate 7 specifying the state of the system. Thus, one of its inherent advantages comes from the fact that 
it is not necessarily restricted to give a vision of the nucleation only focused on the kinetics of clusters. Contrarily, 
a suitable choice of the internal coordinate will lead to the Fokker-Planck equation corresponding to a different 
description of the nucleation process. 

Our purpose in this section will be to particularize this equation for different interpretations of the internal degree 
of freedom. In this way, we will see how we can propose a new Fokker-Planck equation for a global crystallization 
order parameter Q and also how to reproduce the already existing ones. 



A. 7 represents an order parameter. 

A possible way of characterizing the nucleation process, which has recently been used in computer simulations (see 
ref. (3(| , |3?J ) , is by means of an order parameter Q that describes the global degree of crystallization of the system. 
In these simulations, the nucleation rate is calculated from a phcnomcnological rate equation using linear response 
theory. According to our theory, the underlying Fokker-Planck equation describing the kinetics of this process in 
terms of the variable Q, considered as the internal degree of freedom, is given by 

^-M' KQ ' t) ^ + M9 ' l> ') (19) 

where, from eqs. (|l3|) the expressions of the drift and diffusion coefficients are 
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(20) 



If the internal coordinate represents an order parameter associated to the whole crystalinity, the interpretation of 
p is slightly different. Now Q is a variable that characterizes the global state of the system. Therefore, the internal 
space associated to this variable would be the set of all the feasible replicas of the system in a given " macrostate " . 
From this point of view, p would represent the non-equilibrium distribution of this ensemble whose meaning would 
be the fraction of systems with value of the internal coordinate equal to Q at time t. The quantity pdQ would then 
be interpreted as the probability that the order parameter of our system (a single realization of this non-equilibrium 
ensemble) has a value in the interval (Q, Q + dQ) at time t. Since in this case, we work with system replicas that 
obviously do not interact with each other, the system is ideal and eq. (Q) for the chemical potential holds. 



B. 7 related with the number of molecules in a cluster. 

Let us assume that we choose as internal coordinate the number of monomers that constitute a cluster, n, or in 
a more convenient way, the dimcnsionless quantity x — n/N , with N being the total number of monomers. In this 
case p would represent the ( time-dependent ) clusters size distribution, and &(x) would be the free energy needed to 
form a cluster of size x. Assuming that x varies in a continuous way, we would recover the Fokkcr-Planck equation 
(^|) obtained as continuous extrapolation to a master equation, with drift and diffusion coefficients formally analogous 
to expressions ( p0| ) and (pl|). In this case, equation (0) for the chemical potential will be valid only if clusters do not 
interact. This lack of interaction is also implicitly assumed in previous treatments. 



C. 7 is the cluster radius. 

If we consider the coordinate r, corresponding to the position of a monomer in a cluster of radius R, as the internal 
degree of freedom and we work in spherical coordinates, we immediately recover the Smoluchowski equation proposed 
by Ruckenstein and co-workers |2?j 

dp(r,t) _ ( D{f t) y Fp ^ t) + E^A p ^ t )V F G(r)\ (22) 



Notice that now, since the phenomenological coefficient defined in eq. ( |To| ) is in general a function of the internal 
coordinate, the diffusion coefficient will depend on the radial distance r. In this case p would represent the number 
density of monomers at r, and (f> the energy barrier to overcome in order to leave the cluster. 



D. 7 is the set of all system's degrees of freedom. 

If instead of using a single coordinate we carry out a complete description in terms of all the system's degrees of 
freedom Jjj}, characterizing a given configuration, we recover the form of the Fokker-Planck equation proposed by 
Langer f| - @ 

where now D t j — I n this case, p({f]},t) corresponds to the probability density associated to the {r/} 

configuration and $({77}) is the functional free energy. 
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V. THE QUASI-STATIONARY CASE 



In the previous section we have seen how to reproduce some different forms of the Fokker-Planck equations proposed 
in the literature. Our task now will be to ascertain if the equations we have obtained, not only have the appropriate 
form, but are also physically consistent. To this purpose, we will focus our attention on the clusters size distribution 
equation. 

Most of the proposed Fokker-Planck equations are just different continuous diffusive approximations to a discrete 
master equation. Thus, the Tightness of this kind of equations is evaluated in terms of the accuracy of the approx- 
imation. In particular, the two main requirements that a Fokker-Planck equation must satisfy is that it reproduces 
the equilibrium distribution and the deterministic growth rate h plf ; in other words, it must describe nucleation and 
growth appropriately. 

The Fokker-Planck equation we present has been derived directly and not as a continuous approach to a master 
equation. Therefore, in general we do not have an underlying master equation model to which our expression should 
be adjusted and that serve us as a criterion to judge its accuracy. 

However, we have shown that by construction our equation always guarantees the correct equilibrium distribution. 
Moreover, a case exists in which we can build a master equation to evaluate the correctness of our results. It is the 
quasi-stationary case, corresponding to high nucleation barriers. 

When the height of the nucleation barrier is large enough as compared with thermal energy, the system achieves a 
quasi-stationary state characterized by the current 

J(n, t) = J{t) {9{n - ni) - 6{n - n 2 )} (24) 

This form implies that equilibrium is reached independently at each barrier side, consequently the chemical potential 
will be uniform 

n(n,t) = fj.(nx,t)9(no — n) + n(n,2,t)9(n — no) (25) 
Here no specifies the size of the critical cluster. 



Substituting these two equations in (10) and integrating with respect to the coordinate n, we obtain 



D(n ) ( ( g(ng)-j(ni) \ / <S>(n ) - <S>(n 2 ) \\ 
J ^ = I p(ni,t)exp I — I -p(n 2> t)exp\ — I I (26) 

If the nucleation barrier is high enough in such a way that a quasi-stationary state is reached, the clusters distribution 
of intermediate sizes does not vary in time, so that only the states corresponding to p(ni,t) and p(n 2 ,t) matter. The 
value rii represents the smaller cluster distinguishable from equilibrium fluctuations in the metastable phase and n 2 
corresponds to a postcritical stable cluster. Clusters sizes ri\ and n 2 ( where n\ < uq < n 2 ) are chosen such that 
for n < ni, p(n) = Peqin) and for n > n 2 , p{n) — 0. These boundary conditions are habitual in classical treatments 
jl| and final results do not strongly depend on the values of n\ and n 2 pq] . If the intermediate sizes distribution is 
constant, it implies that the lost of one ri\ cluster supposes, after jumping the energy barrier, the formation of one 
postcritical nucleus. Conversely, the gain of a n\ cluster is due to the disappearance of a postcritical cluster. 

Accordingly, the master equation for this quasi-stationary situation has the following form: 



dp(ni,t) 
dt 



k-p(n 2 ,t) - k + p(ni,t) = -J(t) (27) 
dp(n 7 1) 



dt 

dp(n 2 ,t) 
Jt ~ 



0, ni < n < n 2 (28) 



k + p(m,t) - k-p(n 2 ,t) = J{t) (29) 



where k + and fc~ are the forward and backward rates. The forward rate k + is simply the probability that a cluster 
ri\ disappears, or equivalently, that it surpasses the nucleation barrier. Therefore, k + will be proportional to the 
Boltzmann factor associated to the height of the barrier 

fc+ = Aexp(- *<"°>; r *<"-> ) ,30) 
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Similarly, for the rate k one has 



, $(n )-$(n 2 )\ ..... 

Aexp ( kTf ) (31) 



Notice that the constant A enters both expressions of k + and k~ in order to guarantee that in equilibrium the flow 
J eq vanishes. This is the requirement to fulfill detailed balance. 

Therefore, the nucleation rate J(t) obtained from the master equation is given by 

J(t) = A I p(m,t)exp I — I - p(n 2l t)exp I — I I (32) 

Identifying A with the factor in cq. (|2^), we see that at least in the case of high nucleation barriers, our 

diffusive equation reproduces the steady-state nucleation rate. Moreover, we are going to prove that our formalism 
also recovers all the distribution moments of the master equation model. 

By introducing (^) into the continuity equation (|J) we obtain the following expression 

= - J(t) {S(n - m) - 8{n - n 2 )} (33) 

which yields the evolution equation for the r-moment (r = 1, 2, 3 . . .) 

d<n r > 



dt 



= K - rf x )J{t) (34) 



It is immediate to realize that this expression agrees with the corresponding one evaluated from the set of master 
equations d27|)-(f29|). 



Therefore, we have proved the validity of our equation in the quasi-stationary case, in the sense that it satisfies the 
two rightness criteria required to make our approach consistent, and moreover reproduces the same results obtained 
from a master equation. If the height of the barrier is low, the quasi-stationary hypothesis is not longer valid. This 
feature invalidate the considerations we have made in this section. 



VI. CONCLUSIONS 



In this paper we have reviewed the different Fokker-Planck equations proposed in the literature to treat the problem 
of homogeneous nucleation. We have seen that the equations obtained as continuous approximation to a master 
equation have the main drawback of not to be able to reproduce some of the master equation characteristics (as 
the equilibrium distribution, the deterministic growth, the variance of the distribution or the short time behavior). 
Moreover, in all of them (except in the vicinities of the critical point ) the coefficients are not fully determined, but 
rather their expression are postulated, or they are left in terms of the unspecified rates a and /3. To estimate these 
rates an appropriate microscopic kinetic model is required. 

On the other hand, equations obtained from kinetic models base their validity on the correction of the microscopic 
model proposed for the dynamics of the cluster interface. They do not give a complete and correct explicit expression 
for drift and diffusion coefficients. Finally, the field-theory approach is restricted to the vicinity of the critical point 
and to situations where temporal dependencies do not exist. 

We have proposed a new method based upon non-equilibrium thermodynamics in the space of an internal coordinate 
characterizing the state of the system. This method allows us to easily obtain a set of general kinetic equations of the 
Fokker-Planck type that not only reproduce the ones which have been proposed, centered on the evolution of clusters 
size distribution, but they allow alternative descriptions of the nucleation process. In this context, we have derived 
a new Fokker-Planck equation in terms of the crystallinity parameter Q used by Frenkel et al. j36|, (37]] in recent 
simulations. 

We have verified that the Fokker-Planck equation we have obtained apart from having the appropriate form, it 
fulfills the two validity criteria required for this type of equations, at least in the quasi-stationary case. In fact, 
the quasi-stationary case corresponds to low supersaturations ( and slow rates ), that are the typical experimentally 
affordable conditions. 

Our expression also keeps in mind all the possible dependencies (as much in n as in i), and it is able to propose a 
microscopic expression to calculate ( by means simulations ) the drift and diffusion coefficients. 
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Our model allows to treat time-dependent nucleation and also effects of a preexistent clusters distribution. Moreover, 
this new proposed formalism, not only reproduces and overcomes some of the main drawbacks of previous homogeneous 
nucleation treatments, but rather it constitutes an appropriate framework to treat nucleation processes from a wider 
perspective. In this context, we could take into account the possible interactions among clusters just by introducing 
an activity factor in the expression of the chemical potential (fjj). Moreover, the formalism could easily be extended 
to the case of heterogeneous nucleation and even to analyze hydrodynamic effects in the nucleation process by simply 
introducing additional variables in the Gibbs equation characterizing the liquid phase || , || . 
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